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NATIONAL FOREWORD 

This Indian Standard, (Part 3) which is identical with ISO 12213-3: 1997 'Natural gas— Calculation of 
compression factor — Part 3 : Calculation using physical properties' issued by the International Organization 
for Standardization (ISO) was adopted by the Bureau of Indian Standards on the recommendation of the 
Natural Gas Sectional Committee and approval of the Petroleum, Coal and Related Products Division 
Council. 

The text of International Standard has been approved as suitable for publication as Indian Standard 
without deviations. Certain conventions are, however, not identical to those used in Indian Standards. 
Attention is particularly drawn to the following: 

a) Wherever the words 'International Standard' appear referring to this standard, they should 
be read as 'Indian Standard'. 

b) Comma (,) has been used as a decimal marker while in Indian Standards, the current practice 
is to use a point (.) as the decimal marker. 

In this adopted standard reference appears to the following International Standards for some of which 
Indian Standards also exist. The corresponding Indian Standards, which are to be substituted in their 
places, are listed below along with their degree of equivalence for the editions indicated. 



International Standard 

ISO 31-3:1992 Quantities and 
units — Part 3 : Mechanics 

ISO 31-4:1992 Quantities and 
units — Part 4 : Heat 

ISO 6976:1995 Natural gas- 
Calculation of calorific values, 
density, relative density and 
Wobbe index from composition 

ISO 12213-1:1997 Natural gas — 
Calculation of compression 
factor — Part 1 : Introduction and 
guidelines 



Corresponding Indian Standard 



Nil 



Degree of 
Equivalence 



Nil 



IS 1 4504:1 998/ISO 6976:1995 Natural 
gas — Calculation of calorific values, 
density, relative density and Wobbe index 
from composition 

IS 15305 (Part 1):2003/ISO 12213-1:1997 
Natural gas — Calculation of compression 
factor : Part 1 Introduction and guidelines 



Identical 



do 



The Technical Committee responsible for the preparation of this standard has reviewed the provision of 
the ISO 31 -3:1 992 and ISO 31 -4:1 992 and decided that it is acceptable for use in conjunction with this 
standard. 
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Indian Standard 

NATURAL GAS — CALCULATION OF 

COMPRESSION FACTOR 

PART 3 CALCULATION USING PHYSICAL PROPERTIES 



1 Scope 

This International Standard specifies methods for the calculation of connpression factors of natural gases, natural 
gases containing a synthetic admixture and similar mixtures at conditions under which the mixture can exist only as 
a gas. 

This part of ISO 12213 specifies a method for the calculation of compression factors when the superior calorific 
value, relative density and carbon dioxide content are known, together with the relevant pressures and 
temperatures. If hydrogen is present, as is often the case for gases with a synthetic admixture, the hydrogen 
content also needs to be known. 

NOTE — In principle, it is possible to calculate the compression factor when any three of the parameters superior calorific 
value, relative density, carbon dioxide content (the usual three) and nitrogen content are known, but subsets including nitrogen 
content are not recommended. 

The method is primarily applicable to pipeline quality gases within the ranges of pressure p and temperature T at 
which transmission and distribution operations normally take place, with an uncertainty of about ± 0,1 %. For wider- 
ranging applications the uncertainty of the results increases (see annex E). 

More detail concerning the scope and field of application of the method is given in part 1 of this Internationa! 
Standard. 



2 Normative references 

The following standards contain provisions which, through reference in this text, constitute provisions of this part of 
ISO 12213. At the time of publication, the editions indicated were valid. All standards are subject to revision, and 
parties to agreements based on this part of ISO 12213 are encouraged to investigate the possibility of applying the 
most recent editions of the standards indicated below. Members of iEC and ISO maintain registers of currently valid 
International Standards. 

ISO 31 -3:1 992, Quantities and units — Part 3: Mechanics. 

ISO 31 -4:1 992, Quantities and units — Part 4: Heat. 

ISO 6976:1995, Natural gas — Calculation of calorific values, density, relative density and Wobbe index from 
composition. 

ISO 1 221 3-1 : 1 997, Natural gas — Calculation of compression factor — Part 1: Introduction and guidelines. 
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3 Definitions 

Aii definitions relevant to the use of this part of ISO 1221 3 are given in part 1 . 



4 Method of calculation 

4.1 Principle 

The method recommended uses equations which are based on the concept that pipeline quality natural gas may be 
uniquely characterized for calculation of its volumetric properties by an appropriate and distinctive set of measurable 
physical properties. These characteristics, together with the pressure and temperature, are used as input data for 
the method. 

The method uses the following physical properties: superior calorific value, relative density and carbon dioxide 
content. The method is particularly useful in the common situation where a complete molar composition is not 
available, but may also be preferred for its relative simplicity. For gases with a synthetic admixture, the hydrogen 
content needs to be known. 

4.2 The SGERG-88 equation 

The calculation method using physical properties is based on the standard GERG 88 (SGERG-88) virial equation 
for natural gases [^l. 12]. [3]. The standard GERG 88 virial equation is derived from the master GERG 88 
(MGERG-88) virial equation, which is a method of calculation based on a molar-composition analysis W. 

The SGERG-88 virial equation from which the compression factor Z is calculated may be written as 

z = ^ + Bp^ + cp^ ...{1) 

where 

B and C are functions of the input data comprising the superior calorific value Hq, the relative density d, the 
contents of both inert and combustible non-hydrocarbon components of the gas mixture (CO2 and H2) 
and the temperature T; 

pm is the molar density given by 

p^ - p/iZRT) ... (2) 



where 



Z=h{p, T, Hs, d, i:co2' -^Ha) ■ ■ ■ (3) 



However, the SGERG-88 method treats the natural-gas mixture internally as a five-component mixture consisting of 
an equivalent hydrocarbon gas (with the same thermodynamic properties as the sum of the hydrocarbons present), 
nitrogen, carbon dioxide, hydrogen and carbon monoxide. To characterize the thermodynamic properties of the 
hydrocarbon gas adequately, the hydrocarbon heating value Hch is also needed. Therefore, the calculation of Z 

uses 

^ =/2 (p. 7". •'^CH. -^CH. -«:N2. JfC02' -^H2. -*^C0) • • • (4) 

In order to be able to model coke oven gas mixtures, the mole fraction of carbon monoxide is taken to have a fixed 
relation to the hydrogen content. If hydrogen is not present (jchj < 0,001), then set x\-\2 = 0. The natural-gas mixture 
is then treated in the calculation method as a three-component mixture (see annex B). 

The calculation is performed in three steps. 

First, the five-component composition from which both the known superior calorific value and the known relative 
density can be calculated satisfactorily may be found from the input data by an iterative procedure described in 
detail in annex B. 
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Secondly, once this composition is known, B and C may be found using relationships also given in annex B. 

In the third step, equations (1) and (2) are solved simultaneously forpm and Z by a suitable numerical method. 

A flow diagram of the procedure for calculating Z from the input data is shown in figure B.I . 

4.3 Input variables 

4.3.1 Preferred input data set 

The input variables required for use with the SGERG-88 equation are the absolute pressure, temperature and 
superior calorific value (volumetric basis), the relative density, the carbon dioxide content and the hydrogen content. 
Thus the physical properties used the input data set (set A) are 

Relative density is referred to normal conditions (101,325 kPa and °C) and superior calorific value is referred to 
normal conditions (101,325 kPa and °C) and a combustion temperature of 25 °C. 

4.3.2 Alternative input data sets 

Three alternatives to the preferred input data set (see 4.3.1) may be used with the standard GERG virial equation: 

.vNj, Hs, d and xh2 (set B) 
xu2< •tC02' d and ;cH2 (set C) 
xu2< -^002. ■^s and XH2 (set D) 

The alternative input data sets are considered fully in GERG Technical [Ponograph TM5 PI. Use of the alternative 
input data sets gives results which may differ at the fourth decimal place. This part of ISO 12213 recommends the 
use of input data set A. 

4.4 Ranges of application 
4.4.1 Pipeline quality gas 

The ranges of application for pipeline quality gas are as defined below: 



absolute pressure 


OMPa 


^P 


s= 12MPa 


temperature 


263 K 


^ T 


s;338K 


mole fraction of carbon dioxide 





^ XCQ,2 


^0,20 


mole fraction of hydrogen 





^ -*H2 


«0,10 


superior calorific value 


30 MJ m^ 


-3€//s 


=£ 45 MJm-3 


relative density 


0,55 


^d 


s:0,80 



The mole fractions of other natural-gas components are not required as input. These mole fractions shall, however, 
lie within the following ranges: 

methane 0,7 «; xch4 ^1,0 

nitrogen s; x^j ^ 0,20 

ethane *s xcjHg ^0,10 

propane =£= xq^Hs. *= 0,035 

butanes =s XC4H10 ^ 0'0''5 
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pentanes 





^ J^C5H12 


« 0,005 


hexanes 





"^xce 


^ 0,001 


heptanes 





^^Cy 


^ 0,000 5 


octanes plus higher hydrocarbons 





^•«C8+ 


:^ 0,000 5 


carbon monoxide 





^•«C0 


«0,03 


helium 





^XHe 


=£ 0,005 


water 





^^H20 


=£0,000 15 



The method applies only to mixtures in the single-phase gaseous state (above the dew point) at the conditions of 
temperature and pressure of interest. For pipeline quality, the method is applicable over wider ranges of 
temperature and pressure but with increased uncertainty (see figure 1). In the computer implementation, the lower 
temperature limit is set at 250 K. 

4.4.2 Wider ranges of application 

The ranges of application tested beyond the limits given in 4.4.1 are: 



absolute pressure 


OMPa 


^P 


s;12MPa 


temperature 


263 K 


« T 


«338K 


mole fraction of carbon dioxide 





^-^COg 


^0,30 


mole fraction of hydrogen 





« J^H2 


^0,10 


superior calorific value 


20 MJm- 


-3 «//s 


=s 48 MJm-3 


relative density 


0,55 


^d 


<0,90 



The allowable mole fractions of other major natural-gas components are extended to: 

methane 0.5 ^ a:ch4 ^1,0 

nitrogen ^ xj^^ ^ 0,50 

ethane ^ xcjHe ^ 0-20 

propane =£ JccaHe ^ 0,05 

The limits for other minor natural-gas components remain as given in 4.4.1 for pipeline quality gas. 

The method is not applicable outside these ranges; the computer implementation described in annex B will not allow 
violation of the limits of composition quoted here. 

4.5 Uncertainty 

4.5.1 Uncertainty for pipeline quality gas 

The uncertainty in the prediction of the compression factor AZ(for the temperature range 263 K to 338 K) is ± 0,1 % 
at pressures up to lOMPa and ±0,2% between 10MPa and 12 MPa for natural gases with xnj ^ 0,20, 
^C02 =^ 0,09, .rc2H6^0,10 and xh2«0,10, and for 30 MJ m-3 « //s < 45 MJ m-3 and 0,55^^^0,80, (see 
figure 1). 

For gases with a CO2 content exceeding 0,09, the uncertainty of ±0,1 % is maintained for pressures up to 6 MPa 
and for temperatures between 263 K and 338 K. This uncertainty level is determined by comparison with the GERG 
databank on measurements of the compression factor for natural gases [5]. [6] and with the Gas Research Institute 
data [9]. 
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Figure 1 — Uncertainty limits for the calculation of compression factors (Tlie uncertainty limits given are 
expected to be valid for natural gases and similar gases with xf^^ ^ 0,20; jrco2 ^ 0.09; ^C2H6 ^ 0,10 and x^^ ^ 0,^0, 

and for 30 MJ m-3 ?= //g s; 45 MJ m-3 and 0,55 s t/ s 0,80) 

4.5.2 Uncertainty for wider ranges of application 

The estimated uncertainties involved in calculations of compression factors beyond the limits of quality given in 
4.5.1 are discussed in annex E, 

4.5.3 Impact of uncertainties of input variables 

Listed in table 1 are typical values for the uncertainties of the relevant input variables. These values may be 
achieved under optimum operating conditions. 

As a general guideline only, an error propagation analysis using the above uncertainties in the input variables 
produces an additional uncertainty of about ± 0,1 % in the result at 6 MPa and within the temperature range 263 K 
to 338 K. Above 6 MPa, the additional uncertainties are greater and increase roughly in direct proportion to the 
pressure (see reference [3]). 



4.5.4 Reporting of results 

Results for the compression factor shall be reported to four places of decimals, together with the pressure and 
temperature values and the calculation method used (ISO 12213-3, SGERG 88 equation). For verification of 
calculation procedures, it is useful to carry extra digits. 
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Table 1 — Uncertainties of input variables 



Input variable 


Absolute uncertainty 


Absolute pressure 


±0,02MPa 


Temperature 


+ 0,15K 


Mole fraction of carbon dioxide 


+ 0,002 


Mole fraction of hydrogen 


±0,005 


Relative density 


±0,001 3 


Superior calorific value 


±0,06MJm-3 



5 Suppliers of computer programmes 

it is planned to make software available which implements this International Standard. Users are invited to contact 
their ISO member body or ISO Central Secretariat to enquire about the availability of such software. 
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Annex A 

(normative) 

Symbols and units 

The symbols specified in this annex are those which are used in both the main text and in annex B. The units 
specified here are those which give consistency with the values of the coefficients given in annex B. 



Symbol 

/'/YO 
''//I 

%0(1) 

b,,{2) 
R 



CHO 

i^H2 
CHoiO) 

'■ho(2) 
c/yi(2) 

CW2(1) 

^7/2(2) 
c,jki2) 

c 



Meaning 

Zero-order (constant) term in the molar heating value (Wch) expansion of Bn 
[equation (B.20)] 

First-order (linear) term in the molar heating value (Wch) expansion of B-n 
[equation (B.20)] 

Second-order (quadratic) term in the molar heating value (Wch) expansion 
of fiii [equation (B.20)] 

Terms in the temperature expansion of h^Q [equation (B.21)] 
Terms in the temperature expansion of fc//i [equation (B.21)] 
Terms in the temperature expansion of bf^2 [equation (B.21)] 

Terms in the temperature expansion of bjj [equation (B.22)] 

Second virial coefficient [equation (1)] 

Second virial coefficient for binary interaction between component i and 
component; [equation (B.22)] 

Zero-order (constant) term in the molar heating value {Hch) expansion of Cm 
[equation (B.29)] 

First-order (linear) term in the molar heating value (Hch) expansion of Cm 
[equation (B.29)] 

Second-order (quadratic) term in the molar heating value (Wch) expansion 
of Ciii [equation (B.29)] 

Terms in the temperature expansion of cfio [equation (B.30)] 
Terms in the temperature expansion of c//i [equation (B.30)] 
Terms in the temperature expansion of ch2 [equation (B.30)] 

Terms in the temperature expansion of eg [equation (B.31)] 
Third virial coefficient [equation (1)] 



Units 

mS-kmoM 

m3.MJ-i 

m3kmol MJ-2 

m3 kmol-i 
m3kmol-i K-1 
m3kmol-i K~2 

m3.MJ-i 

m3.MJ-i.K-i 

m3.MJ-iK-2 

m3kmolMJ-2 
m3.kmolMJ-2 K-1 
mS-kmoi MJ-2 K-2 

m3kmol-i 
m3.kmol-i K-1 
m3.kmol-i K-2 

m3.kmol-i 
m3-kmol-i 

m6kmol-2 

m^ kmoH MJ-1 

m6.MJ-2 

m6kmol-2 

m6.kmol-2K-i 

m6kmol~2K-2 

m6kmol-iMJ-i 
m6.kmol-i MJ-1 K-1 
m6.kmol-iMJ-iK-2 

m6.MJ-2 

m6MJ-2.K-i 

m6.MJ-2.K-2 

m6kmol-2 

m6.kmol-2.K-i 

m6kmol-2.K-2 

m6kmol-2 



IS 15305 (Part 3) : 2003 
ISO 12213-3 : 1997 



Symbol Meaning 



d 
// 

M 

P 

R 

T 

t 

V'm 

X 

V 



Z 
P 

Pn 



Third virial coefficient for ternary interaction between components ij and k 
[equation (B.31)] 

Relative density [d{a\r) = 1 ; equation (B.1)] 

Change in the molar heating value Hqh dunng iteration 
[equations (B.10) and (B.11)] 

Superior calorific value [gas at normal conditions (0 °C, 1,013 25 bar), 
combustion temperature 25 °C] 

Molar heating value (combustion temperature 25 °C) 

Molar mass [equations (B.5) and (B.8)] 

Absolute pressure 

(Universal) gas constant 

Absolute temperature 

Celsius temperature [= r- 273,15; equation (B.27)] 

Molar volume (= 1/pm) 

Mole fraction of a component 

Combination rule parameters for the binary unlike-interaction virial coefficients 
Bt2 and 5-13 (table B.I) and the ternary unlike-interaction virial coefficient C,;^ 
[equation (B.32)] 

Compression factor 

Mass density [equations (B.8) and (B.42)] 

Molar density (= V^-'') 



Units 

m^kmohs 

MJkmol-1 
MJ.m-3 

MJkmol-1 

kgkmol-"' 

bar 

m3barkmohiK-i 

K 

°C 

m3kmol-i 



kg-m-3 
kmolm-3 



Additional subscripts 

n Value at normal conditions (^n = 273,15 K, p^ = 1,013 25 bar) 

CH For the equivalent hydrocarbon 

CO For carbon monoxide 

CO2 For carbon dioxide 

H2 For hydrogen 

N2 For nitrogen 



Additional qualifiers 

(air) For dry air of standard composition [equation (B.I)] 

(D) For special value of p used in equation (B.1 1 ) 

1 For the equivalent hydrocarbon [equations (B.I 2) and (B.I 5)] 

2 For nitrogen [equations (B.I 2) and (B.I 6)] 

3 For carbon dioxide [equations (B.12) and (B.17)] 

4 For hydrogen [equations (B.12) and (B.I 8)] 

5 For carbon monoxide [equations (B. 1 2) and (B.1 9)] 
(id) Ideal gas state 

(u) Iteration counter (B. 2.1) 

(i') Iteration counter (B.2.2) 

(h) Iteration counter (B.4) 
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Annex B 

(normative) 

Description of the SGERG-88 method 

This annex gives the equations for and numerical values of coefficients which together specify completely the 
SGERG method for calculation of compression factors. 

It also deschbes iteration procedures adopted by GERG [3] for implementing the method in the verified Fortran 77 
subroutine SGERG. FOR. This subroutine provides the correct solution; other computational procedures are 
acceptable provided that they can be demonstrated to yield identical numerical results. The calculated results shall 
agree to at least the fourth place of decimals with the examples given in annex 0. 

Other implementations which are known to produce identical results are as follows: 

a) A Basic version, described in GERG TM5 Ol, which may be used with a variety of metric reference conditions. 
This programme was designed mainly for PC applications. 

b) A version in C, described in German DVGW Directives, sheet G486 [8]. 

c) A version in Turbo Pascal. 

Ah these programmes have Been verified to give the same results to within 10-5. The availability of the programmes 
and the conditions which apply to their use are discussed in part 1 of this International Standard. 



B.1 Basic structure of the calculation method 

As described in 4.2, the calculation proceeds in three steps, which are shown schematically in figure B.1. 



Input data: 



p. t, W5, d, x^oj, x„2 



Calculate 


intermediate 


data; 




T„ -273,1SK 






g„.x^Q, X 


h(u), X 


u-Au). 


H,hM 






Mi„(u). eM 


B.M 


Pm.r 


M. HsM 
(see clause B.2) 








I 













Calculate virial coefficients; ?" = (/♦ 273,15) K 
8(1 Hr„, xX CiT. WfH. ^/) 

(see clause B,3) 
















Solve virial equation: 

Z-^*Bp^^T,p. x,U CeJiT.p. x,) 

(see clause B. 4) 





Figure B.I — Flow diagram for standard GERG'88 calculation method 

{xi = mole fraction of component i) 
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The calculation is described below in the order in which these three steps are carried out. 

Step I 

The input data are pressure, temperature, gross calorific value, relative density and the mole fractions of carbon 
dioxide and hydrogen. If the values of the first three parameters are in any units other than bar, °C and MJ/m3, they 
shall first be converted precisely to values in bar, °C and MJ/m^, respectively, using the guidelines set out in 
annex D. 

The input data are then used to calculate the following intermediate data: 

Mole fraction of: 

hydrocarbon gas jcch 

nitrogen xf^^ 

carbon monoxide j^co 

Molar heating value of the equivalent hydrocarbon //ch 

Molar mass of the equivalent hydrocarbon AfCH 

Second virial coefficient (Tn = 273,1 5 K) Bn 

Molar density at normal conditions Pm,n 

Mass density at normal conditions Pn 

Superior calorific value of the gas //s 

In equations (B.I) to (B.46), each symbol represents a physical quantity divided by its selected unit (see annex A), 
such that their quotient is the dimensionless value of the quantity. 

Step II 

The intermediate data are used to calculate the second and third virial coefficients for the natural gas at the required 
temperature, B(7,//cH.^.) a"d C{T,HcH.^,)- 

Step 111 

The second and third virial coefficients determined in the second step are inserted in the virial equation, and the 
compression factor Z is calculated for a given pressure and temperature. 

The symbols used are defined in annex A. 



B.2 Calculation of intermediate data 

The eight intermediate-data values (xch- -^N2' -^co. ^ch. ^ch. ^n. Pm,n. Pn) are determined from equations (B.1) to 
(B-8) using the iterative method presented in figure B.2, Values of the constants used in these equations are given 
in table B.1 . 

Pn-^Pn(air) ...(B.1) 

^co = 0,096 4 JCH2 • ■ ■ (B.2) 

Vn.,n(id)=RTjp, ...(B.3) 

An,n(v) = [Vn(W)+fin(v)]"' • • • (B.4) 



10 



input data: 



A /. Ws.rf, Xro2. ■><■« 



Calculate; 



Pn.JfcO 



Equations (B.1), (B.2) 



u = 0, 1/= 







t 




«CH 


((; = 


0) = 


1000 


Sr 


(v 


= 0) = 


- 0.065 



Calculate: 



em,n^ = o) 



Equations <B. 3). (B. 4; 



u = u ■>■ 1 

1 



Calculate: 


X(h{u), Xfijiu). M^Ju). Pniu) 

Equations (B.5) to (B.8) 



Calculate: 


r„ = 273,15 K 
Equations (B. 12) to (B. 25) 



No 




Result: 



•"■CH. •*'n2' -*'cO' '^CH 
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W[h(c/) = H[h'"- "" * DH^„l.u) 
Equation (B.10) 



Calculate: 




Figure B.2 — Flow diagram for computing intermediate data by iteration 



11 
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Mcni'i) = - 2,709 328 + 0,021 062 199 «ch(« - 1) 

-*N2 (") = 1 - -^Ch(") - ->^C02 - -^Hs - -^CO 

Pn(«) = {xCHiu)McH(") + -^N2 («)WN2 ] Pm,nM + (vc02^C02 + ^H2'^H2 + '^CO'^Cojpm.ni^ 



. . (B.5) 
■ ■ (B.6) 
• ■ (B.7) 
. . {B.8) 



Table B.1 — Values of the constants used in equation (B.I) 

(adjusted to conform with the molar masses and molar calorific values in ISO 6976:1975) 



//H2 


= 285,83 


MJkmol-1 


//CO 


= 282,98 


MJkmol-1 


Nf^2 


= 28,013 5 


kg kmol-i 


^C02 


= 44,010 


kg-kmol-l 


Mhz 


= 2,015 9 


kgkmol-i 


Mqo 


= 28,010 


kgkmohi 


R 


= 0,083 145 1 


m3.barkmol-l K-1 


Vm,n('^) 


= 22,414 097 


m3 kmoj-l 


Pn(air) 


= 1,292 923 


kg m-3 



B.2.1 Iteration with the molar heating value Hc\-\ (inner loop) 

Equations (B.I) to (B.8) are applied in sequence so as to obtain the first approximation in the wthr iteration step. The 
starting values are: 

//ch(" = 0) = 1 000 MJkmol-1 

Br,{v == 0) = - 0,065 m3.kmoh1 

The values of the other constants used in equations (B.I ) to (B.8) are given in table B.I . 

The convergence criterion for this inner iteration loop is that the absolute difference between the calculated density 
of the gas at normal conditions pn(u) and the known density (either measured directly or determined from the 
relative density) of the gas at normal conditions pn is less than 10-6, j.e. 



|Pn -Pn(«)|<10~^ 



. . (B.9) 



If this condition is not satisfied, then an improved value of the molar heating value Hch(li), for use in equations (B.5) 
to (B.8), is calculated using equation (B.10) as follows: 



//ch(«} = /^ch(«-1)+£»//ch(«) 
where 

/Wch(") = [Pn - Pn(«)] [piD) - Pn(«)] 



-1 



(B.10) 



.(B,11) 



pr,(u) being the density value for the current iteration step [commencing with //ch(w - "•)]. 

p{D) being the density determined by equations (B.4) to (B.8) using [//ch(« - 1) + 1] as input for the molar 
heating value. 
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When the left-hand side of equation (B.9) is less than 10-6, this iteration loop is terminated and iteration with the 
second virial coefficient begins. 

B.2.2 Iteration with the second virial coefficient Bn (outer loop) 

The intermediate values xcni"), ■rN2("). Jcco and //ch(") ffom *he preceding iteration and the input data j:co2 and jthj 
are used to determine an improved value for the second virial coefficient Bn(v) for the whole gas at normal 
conditions. 

The second virial coefficient for the natural gas is given by the following equation: 

S(T) = x?Bii +2x1^2^12 +2x-^X3B^3+2x-^X4B^^ +2x^x58^5 +4^22 +2x2x3823+2x2X4624 + 

+ 4B33 + x2fi44 + xffiss . . . {B.12) 

Some of the terms that are missing in equation (B.12), i.e. B25, 634, etc., have been found not to improve the 
accuracy of the calculation if included and are therefore set at zero. 

Bn(v)=fi(7-n) . . . (B.13) 

where 

r = 77, =273,15 . . . (B.14) 

•«i=J^ch(«) ...(B.15) 

X2=xu2{u) ...(B.16) 

^3=J^C02 ...{B.17) 

^4=^H2 ...{B.18) 

A-5=^C0 ...(B.19) 

Bn = bno + /'//I^ChC") + f>H2"CH^M ■ ■ ■ (B-20) 

where the coefficients bf^Q, b^-\ and bn2 are second-degree polynomials as a function of temperature 

fii 1 - ^«ofO) + bt^oi^)T + b„o(^)T^ + Ki^OJ + I'mOF + bf^,(2)T^]HcH(u) + 

+ [bH2{0) + bf,2m + bH2(2)T^H^y^{u) ...(B.21) 

and the second virial coefficients B14, 615, 622- ^23. ^24. ^33. ^34. B44 and 655 are also second-degree polynomials 
as a function of temperature, in the general form 

Bij = bij (0) + b^j (1)7- + bij (2)r2 ... (B.22) 

The unlike-interaction virial coefficients B-^2 and B13 are given by 

B12 = [0,72 + 1,875x1 0-5 (320- 7)2 ](Bii + B22)/2 . . . (B.23) 

Bl3 = - 0,865(fiiifi33)^^^ - . . (B.24) 
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The coefficients in equations (B.21) to (B.24) are given in table B.2. 

Table B.2 — Numerical values for the coefficients b(0), b{^) and b(2) in the temperature expansion of the 
second virial coefficient for pure gases and of the unlike-interaction virial coefficients 

(Tine units of B are m3kmol-i when the temperature is in kelvins) 







y 


b{0) 


b(\) 


b[2) 




CH 




HO 


-4,254 68x10-1 


-2,865 00x10-3 


-4,620 73x10-6 




CH 




m 


8,771 18x10-4 


-5,562 81 X10-6 


8,815 10x10-9 




CH 




H2 


-8,247 47x10-7 


4,314 36x10-9 


-6,083 19x10-12 




N2 




22 


-1,446 00x10-1 


7,409 10x10-4 


-9,119 50x10-7 




CO2 




33 


-8,683 40x10-1 


4,037 60x10-3 


-5,165 70x10-6 




H2 




AA 


-1,105 96x10-3 


8,133 85x10-5 


-9,872 20x10-8 




CO 




55 


-1,308 20x10-1 


6,025 40x10-4 


-6,443 00x10-7 




CH + 


N2 


12 


y =0,72 + 1,875 X 10-5(320- 


-7)2 






CH + 


CO2 


13 


y = - 0,865 






« 


CH + 


H2 


14 


-5,212 80x10-2 


2,715 70x10-4 


-2,500 00x10-7 




CH + 


CO 


15 


-6,872 90x10 2 


-2,393 81 X 10-6 


5,181 95x10-7 




N2 + 


C02 


23 


-3,396 93x10-1 


1,611 76x10-3 


-2,044 29x10-6 




N2 + 


H2 


24 


1.200 00x10-2 


0,000 00 


0,000 00 





The value of B^{\>) obtained from equation (B.13) is used to calculate the vth approximation of pmn using equation 
(B.4). 

Equation (B.6) is then used, in the inverse way to that in which it was used previously, to obtain a value for Hq{v), i.e. 

Ws(0 = h(")WcH("-1) + JC4/^4 +J^5W5]An,nM . . . (B.25) 

where H^{= H^^ and Hs(=Hco) are the molar heating values at 298,15 K of hydrogen and carbon monoxide, 
respectively. The convergence criterion for the outer iteration loop (iteration counter v) is that the absolute difference 
between the measured superior calorific value Hs and the calculated calorific value Hs(v) is less than 10-4, i.e. 

\Ms-Ms{v)\<^0-^ ...(B.26) 

If this criterion is not satisfied, then the value for BnM determined from equation (B.13) is used as a new input value 
for equation (B.4) and the whole iteration procedure, i.e. the inner iteration loop (iteration counter «), is restarted 
from equation (B.5) using the current values of //ch(" - 1) and Pm.niv)- 

When both convergence criterion (B.9) and convergence criterion (B.26) are satisfied simultaneously, the final 
intermediate data for the mole fractions xch and xn2 and for the molar heating value Wch have been determined. 



B.3 Calculation of virial coefficients 

The second and third virial coefficients 8(7) and C(7) of a natural gas are now determined from the mole fractions 
1CO2 and XH2 (input data) and xch, xh2 and xqq (intermediate data) and the molar heating value Hch (see figures 
B.1 and B.3). 
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Calculate: 



Input data 
P 


(th 

t. 


rd iteration loop); 


i 


w = 


1 


Calculate; 






r = (/*2?3.15K) 


1 


Calculate; 






BiT). an. ejo) 












1 


w - w * ^ 


i 


Calculate; 






Equations (B.38), (B.39) 




Z{p, T) 



Equation (B.41) 



Figure B.3 — Flow diagram for compression factor calculation 



B.3.1 Calculation of 8(7) 

The second virial coefficient B(T) is calculated from equation (B.12) by the procedure described in B.2.2 for a 

temperature 



r = r + 273,15 



. . . {B.27) 



B.3.2 Calculation of 0(7) 

The third virial coefficient for a natural gas at a temperature 7 is determined using the following equation: 

C{T) = xfC-^i^ + 3a-2x2Ci-,2 + 3jff JC3C113 + 3^:^X4^14 + 8x^2^5^15 + 3xijc|Ci22 + 6x1x2^30123 + 

+ 3xix§Ci33 + x|C222 + 3x|x3C223 + SxsxfQaa + xfCaaa + x|C444 . . . (B.28) 

The possible additional terms that are missing in equation {B.28) have been found not to improve the accuracy of 
the calculation if included and have therefore been set at zero. 
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Furthermore, in equation (B.28); 

C"i 1 1 = Cho + LHiWcH + <'H2/^CH ■ ' " ^^'^^^ 

where c^q, chi and ch2 are second-degree polynomials as a function of temperature, viz: 

Cm = rHo(O) + CHo(1)r+ f-Ho (2)7-2 + ^,.^^(0) + .^^{l)?- + c^^{2)T^]Hch + 

+ [rH2 (0) + .■H2 0)T + <H2 (2)^^^ ]« 2^ . . . {B.30) 

as are C222, C333, C444, C115, C223 and C233, viz: 

C,j, = c,, (0) + c,j, (1)7 + c^j, (2)^2 . . . (B.31) 

The coefficients in equations (B.30) and (B.31) are given in table B.3. 

Table B.3 — Numerical values of the coefficients r(0), c(1) and c(2) in the temperature expansion of the third 
viriai coefficient for pure gases and of the unlike-interaction virial coefficients 

(The units of C are m6.kmol-2 when the temperature is in kelvins) 





ijk 


c(0) 


r(1) 


r(2) 


CH 


HO 


--3,024 88 X 10-' 


1,958 61 X 10-3 


-3,16302x10-6 


CH 


H^ 


6,464 22 X 10 -l 


-4,228 76 X 10-6 


6,881 57 X 10-9 


CH 


H2 


-3,328 05 X 10-7 


2,231 60 X 10-9 


-3,677 13 X 10-12 


N2 


222 


7,849 80x10-3 


-3,989 50 X 10-5 


6,118 70x10-8 


CO2 


333 


2,051 30x10-1 


3,488 80x10-5 


- 8,370 30 X 1 0-8 


H2 


444 


1,047 11 x 10^3 


-3,648 87 X 10-6 


4,670 95 X 10-9 


CH + CH + N2 


112 


>• = 0,92 + 0,001 3(7-270) 






CH + CH + CO2 


113 


y = 0,92 






CH + CH + H2 


114 


y = 1 ,20 






CH + CH + CO 


115 


-7,367 48x10-3 


-2,765 78 X 10-5 


3,430 51 X 10-8 


CH + N2 + N2 


122 


y = 0,92 + 0,001 3 (T - 270) 






CH + N2 + CO2 


123 


>•= 1,10 






CH + CO2 + CO2 


133 


y = 0,92 






N2 + N2 + CO2 


223 


5,520 66 X 10-3 


- 1,686 09x10-5 


1,571 69x10-8 


N2 + CO2 + CO2 


233 


3,587 83 X 10-3 


8,066 74 X 10-6 


-3,257 98x10-8 



The other unlike-interaction virial coefficients used are given by 

Cijk = yijk (CiuCjjjC 1,1,1, ) ... (B.32) 
where yiji, is given by 

V1 12 = .V122 = 0,92 + 0,001 3(7 - 270) . . . (B.33) 

>H3=>'133=0,92 ...{B.34) 

vi 14 =1.20 ...(B.35) 

vi23-1,10 ...(B.36) 

Equation (B.32) shows that the temperature dependence of the unlike-interaction virial coefficients is determined 
essentially by the temperature dependence of the third virial coefficients for the pure components. 
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B.4 Calculation of the compression factor and molar density 



The very last stage in the calculation of the compression factor and the molar density is to solve equations (1) and 
(2) simultaneously tor the given value of the pressure p. For the first approximation in the iteration using w, pm is 

given by 

P^\w^O)^RTIp + B ...(B.37) 

where the second virial coefficient B is defined by equation (B.I 2) for a temperature T (see figure B.3). An improved 

value prn('tO is then given by 

p^'(u) - {t<rlp{^ + Bp^iw - 1) + Cp^iw - 1)] ... (B.38) 

where the third virial coefficient C for the mixture is defined by equation (B.28) for a given temperature T. The 
convergence criterion for the iteration using w is that the absolute difference between the calculated pressure p(w) 
given by equation (B.39) and the given pressure /j is less than 10-5 [see equation(B.40)]. 

p{^v) = /?rp^(vv)[l + Bp^{w) + Cp^(w)] . . . (B.39) 

|/)-/7(vv)|<10-5 ...(B.40) 

If this condition is not satisfied, then the current value for the molar density pmiw) is used as the new value pm{w- 1) 
in equation (B.38) and an improved value of the molar density pm(w) is calculated. 

However, if the left-hand side of equation (B.40) is less than 10-5, the iteration routine is ended, and Pm(w') is the 
final molar density pf^. The compression factor is then given by 

Z = 1 + /?/^+Cp2 ...(B.41) 

NOTE — The mass density can be calculated as follows; 

p = [<Wa\r)pZ^ ^n/(/'n2^)] • • • (B-42) 

Z and Zp being rounded to four places of decimals before being used in the density calculation. 
Report the density to 3 significant figures. 



B.5 Consistency checks on the SGERG-88 method 

The following tests, which provide partial consistency checks on the input data, shall be applied when carrying out 
calculations by the SGERG method. 

a) The input data shall satisfy the following condition: 

d > 0,55 + 0,97,vco2 - O.45.XH2 . . . (B.43) 

b) The intermediate calculated value for the mole fraction of nitrogen shall satisfy the following conditions: 

- 0,01 s= x^^ =s 0,5 ... (B.44) 

-^N2 +-^C02 =^0,5 ...(B.45) 

c) Furthermore, the internal consistency of the input data for the third iteration loop shall satisfy the condition: 

d > 0,55 + 0,4x^2 + 0,97jcco2 ~ 0,45ah2 . . . (B.46) 
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Annex C 

(normative) 

Example calculations 

The following example calculations shall be used for the validation of computer implementations of the SGERG 
method not cited in annex B. The calculations have been carried out using the validated executable programme 
GERG88.EXE, which incorporates the subroutine SGERG. FOR described in annex B. 

Table C.I — Input data 





Gasi 


Gas 2 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


->^C02 


0,006 


0,005 


0,015 


0,016 


0.076 


0,011 


-rH2 


0,000 


0,000 


0,000 


0,095 


0,000 


0,000 


d 


0,581 


0,609 


0,650 


0,599 


0,686 


0,644 


Hs (MJ m-3) 


40,66 


40,62 


43,53 


34,16 


36,64 


36,58 



Table C.2 — Results (Z-values) 



Cond 


itions 


Gas 1 


Gas 2 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


bar 


°C 














60 


-3,15 


0.840 84 


0,833 97 


0,794 15 


0,885 69 


0,826 64 


0,854 06 


60 


6,85 


0,862 02 


0,856 15 


0,822 10 


0,901 50 


0,850 17 


0,873 88 


60 


16,85 


0,880 07 


0,875 00 


0,845 53 


0,915 07 


0,870 03 


0,890 71 


60 


36,85 


0,908 81 


0,904 91 


0,882 23 


0.936 84 


0,901 24 


0,917 36 


60 


56,85 


0,929 96 


0,926 90 


0,908 93 


0,953 02 


0,923 94 


0,936 90 


120 


-3,15 


0,721 46 


0,711 40 


0,643 22 


0,808 43 


0,695 57 


0,749 39 


120 


6,85 


0,759 69 


0,750 79 


0,690 62 


0,836 13 


0,738 28 


0,784 73 


120 


16,85 


0,792 57 


0,784 72 


0,731 96 


0,859 99 


0,774 63 


0,814 90 


120 


36,85 


0,844 92 


0,838 77 


0,797 78 


0,898 27 


0,831 66 


0,862 66 


120 


56,85 


0,883 22 


0,878 32 


0,845 54 


0,926 62 


0,872 69 


0,897 49 



These gases are the same as the six gases in ISO 12213-2, annex C, where the complete molar compositions are 
given. 
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Annex D 

(normative) 

Conversion factors 

D.I Reference conditions 

The reference conditions for which the standard GERG-88 virial equation was developed and which the 
SGERG.FOR computer subroutine uses internally are 



Calorific value by combustion at 



7i =298,15K(fi =25°C) 
;>=101,325kPa 



gas metered at ^2 ^^^ 273,15 K (ra = °C) 

p = 101,325 kPa 

The latter set of conditions are also the reference conditions for relative density. 

Considerable care is needed to ensure that correctly referenced inputs are used for calorific value and relative 
density. Several countries normally use the above conditions, but others use alternative conditions. This can easily 
cause confusion, particularly since the unit of calorific-value measurement in each case may still be MJ m-3. 
Table D.I is a guide to which of the major international gas-trading countries use which reference conditions. 

For those using non-metric units for calorific value (i.e. Btu ft-3), conversion both of units and of reference conditions 
is required. 

The conversion factors used are taken from reference [3]. 

Table D.1 — Nationally adopted metric reference conditions for the measurement of calorific value 





'1 rc) 


>2 (°C) 


Australia 


15 


15 


Austria 


25 





Belgium 


25 





Canada 


15 


15 


Denmark 


25 





France 








Germany 


25 





Ireland 


15 


15 


Italy 


25 





Japan 








Netherlands 


25 





Russia 


25 


0or20 


United Kingdom 


15 


15 


United States of America 


15 


15 


NOTES 




1 In all countries the reference pressure is 101,325 kPa {= 1,013 25 


bar). 


2 ?i is the combustion reference temperature. 




3 t2 is the gas-metering reference temperature. 
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D.2 Units and conversion factors for pressure and temperature 

if the input variables p and ; are not in the necessary units of bar and °C, then conversions must be mad© in order to 
use the Fortran implementation. A selection of appropriate conversion factors is given in table D.2. 



Table D.2 — Conversion factors for pressure and temperature 



Pressure 


p (bar) 


= [p(kPa)]/100 


P (bar) 


= b(MPa)]x10 


p (bar) 


= [p(atm)]x 1,013 25 


p (bar) 


= [p(psia)]/14,503 8 


p (bar) 


= b(psig) + 14,695 9]/14,503 8 


Temperature 


r (°C) = 


r(K)- 273,15 


t (°C) = 


[r(°F)-32]/1,8 


' rc) = 


[/(°R)]/1,8-- 273,15 



D.3 Units and conversion of calorific value and density between reference conditions 

Because both superior calorific value and relative density are functions of the composition of a gas mixture, and 
because the thermophysicai properties of the individual components depend upon temperature and pressure in 
individual ways, it is in principle impossible (without knowledge of the composition) to convert the calorific value and 
the relative density, known at one set of reference conditions, to exact corresponding values for any other set of 

reference conditions. 

However, because the relevant reference conditions are always thermodynamicaily close together, and because 
natural gases do not vary in their composition to any major extent, it is possible in practice to give conversion 
factors which may be applied to any typical natural gas with essentially no loss of accuracy in the converted 
physical properties. 

D.3.1 Conversion factors for the units 

If the input variable Hq is not in the necessary units MJ m-3, then a conversion must be made. 

The appropriate conversion factors for calorific values at a given set of reference conditions (see clause D.I) are 
given in table D.3. 



Table D.3 — Conversion factors for calorific value 



Ws (MJ m-3) = [//s(kWh m-3)] x 3,6 
Hs (MJ m-3) = [Hs(Btutt-3)]/26,839 2 
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D.3.2 Conversion between different reference conditions 



If the input variables H^ and d are not at the correct reference conditions (see clause D.1), then a conversion nriLtst 
be made. Table D.4 gives the conversion factors appropriate to the reference conditions in table D.1 . 



Table D.4 — Conversion factors for calorific value and relative density 



Calorific value //g at reference conditions 




t^ = 


25 °C, f2 


= 0°C,p2- 1,013 25 bar: 




Hs 


= Ws(fi = 


0°C,r2= 0°C,p2^ 


1,013 25 bar) X 


0,997 4 


Hs 


= //s(fi = 


15°C,f2 = 15°C,p2 = 


1,013 25 bar) X 


1,054 3 


Hs 


= Hs(h = 


60 °F, f2 = 60 °F, pz = 


1,015 92 bar) X 


1,053 5 


Hs 


= Ws(fi = 


60 °F, r2 = 60 °F, pz = 


1,015 60 bar) X 


1 ,053 9 


Relative density d at reference conditions 




^2 = 


--0°C,p2-- 


= 1,013 25 bar: 






d = 


d(t2 = ^5 


°C,;72 = 1,013 25 bar) 


X 1,000 2 




d^ 


d [t2 = 60 


°F, p2 = 1,015 92 bar or 1 ,01 5 60 bar) 


X 1,000 2 


NOTE — p -- 


= 1,015 60 bar =14,73 


psia (pressure 


customary in US). 
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Annex E 

(informative) 

Performance over wider ranges of application 

The standard GERG virial equation has been comprehensively tested, over the temperature range 263 K to 338 K 
and at pressures up to 12 IVIPa, with the GERG databanl< [5] and the Gas Research institute data [9] for gases within 
the ranges of composition, superior calorific value and relative density given for pipeline quality gases {see 4.4.1). 
Within these limits, the uncertainties are as given in 4.5. 

Rough estimates of the uncertainties involved in calculations of compression factors for wider ranges of application 
are plotted in figures E.1 to E.4 as pressure-composition plots for nitrogen, carbon dioxide, ethane and propane, 
respectively. 

In figures E.I to E.4, the performance of the SGERG equation is illustrated up to a maximum pressure of 30 MPa. 
This is for informative purposes only and is not intended to imply a recommendation that the equation be used 
routinely or uncritically above the normal specified limits. The uncertainty limits are dependent upon pressure, 
temperature and composition, and are also strongly affected by the proximity of the phase boundary. The estimated 
uncertainty limits presented below are based upon less comprehensive data, published as a supplement to the 
GERG databank l^l, and upon the databanks in references [5] and [9]. 






30 



20 



10 











SGERG-8e equation 

7'=263Kto338K 




- 




D 
















^^-----.^ 




A 
B 
C 
D 


AZ -iO.1% 
±0,1% to ±0,2% 
±0,2 % to ±0,5 % 
±0,5% to ±3,0% 

1 1 




^^^^^~--^° 




A 

1 1 





0,1 



0,2 



0,3 



0,U 0,5 

Mole fraction of nitrogen, Xn2 



Figure E.1 — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high nitrogen content 
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The worst-case limits have always been chosen. Dashed lines are used to separate two regions of estimated 
uncertainty when the experimental evidence is not sufficient to determine the position of the boundary. The detailed 
composition of the gas will have a strong influence on the position of the phase boundary and the user should, 
therefore, make his own phase boundary calculation. 



- 30 



20 



10 



S6ERG-8e equation 

r = 263 K to 338 K 




A AZ=!±0,1% 

B ±0,1%roi0,2% 

C ±0,2% to ±0,5% 

D iO,5%tOi3,0% 



0,1 



0,2 0.3 

Mole fraction of carbon dioxide, x^qj 



Figure E.2 — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high carbon dioxide content 
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S 30 



20 



10 



A AZ < ±0,1 % 

B ±0,1 % to tO,2 % 

C i0,2 % to ±0,5 % 

D ±0,5 % to ±3.0 % 



SGERG-88 equation 

r = 263Kto338K 




0,05 



0,1 



0,15 0,2 

Mole fraction of ethane. a'c2H6 



Figure E.3 — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high ethane content 
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i 30 



20 




10 



S&ERG-88 equation 

r = 263Kto338K 



A A/^tO,1% 

B ±0,1% to ±0,2% 

C ±0,2% to ±0,5% 

D ±0,5 % to ±3,0 % 



0.05 



0,1 
Mole fraction of propane, x,; 



Figure E.4 — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high propane content 



The overall results at pressures up to 10MPa and temperatures within the range 263 K to 338 K can be 
summarized as follows. Only gases having mole fractions within the limits given below will have uncertainties within 
± 0,1 %, ± 0,2 % and ± 0,5 %, respectively, within the given pressure and temperature domain. 



Component 


Mole fraction for an uncertainty within | 




±0,1 % 


± 0,2 % 


+ 0,5 % 


Nitroger\ 


«0,20 


«0,50 


— 


Carbon dioxide 


«0,09 


=s0,12 


«0,23 


Ethane 


«0,10 


^0,11 


=^0,12 


Propane 


« 0,035 


^0,04 


^ 0,045 
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Annex F 

(informative) 

Subroutine SGERG.FOR in Fortran 

c 

Q* ************************* -I, ***** -k * -k * -k ■), -k i, * * * i, ***** -k * * -k ********* i, * i, * -k * * * 1, * i, 

C last update: 07.08.96 M. Jaeschke, J. Sikora 

C 

Q ******************* GERG-8 8 VIRIAL EQUATION *********************** 

C SUBROUTINE SGERG 

C 

C 'SGERG' CALCULATES THE COMPRESSION FACTORS OF NATURAL GASES USING 

C A SIMPLIFIED GAS ANALYSIS 

C 

C LEGAL COPIES OF THIS PROGRAM MAY ONLY BE OBTAINED FROM THE 

C MEMBERS OF THE GERG WORKING GROUP ON COMPRESSION 

C FACTORS OF NATURAL GAS AS GIVEN IN THE GERG TECHNICAL 

C MONOGRAPH TM 5 (1991) . 

C 

C The calculations are based on the following four 

C input parameters for the gas analysis: 

C 

C 

c 
c 
c 
c 

C *) note : metering atT = O.OC , P= 1.01325 bar 

C : coinbustion at T = 25.0 C 

C 

C Further input parameters used are: 

C 

C P : pressure in bar ( i-> 120 ) 

C TC: temperature in degrees Celsius (-23 -> 65 ) 

C 

C A CALCULATED value for X2 , the molar fraction of nitrogen 

C is returned. 

C 

C The calculated values are: 

C 

C Z : compression factor 

C D : molar density in mol/m**3 

C 

C **** pqj- gome compilers the SAVE option has to be set explicitly **** 

C 

C SYNTAX : CALL SGERG (X2 , X3 , HS , RM, X5 , P, TC, Z, D) 

C 

C The coefficients used in this program are conform with the 

C values given in subroutine GAS682, from July 20,1988, appended 

C to report 8807, Van der Waals Laboratory, Amsterdam. 

C 



-1- 


X3 


-2- 


HS 


-3- 


RM 


-4- 


X5 



mole fraction C02 
calorific value in MJ/m'^3 *) 
relative density *) 

mole fraction H2 



Valid ranges: 
0.0 -> 0.3 
20 -> 48 

0.55 -> 0.9 
0.0 -> 0.1 



J.P.J. Michels & J. A. Schouten 
August 16, 1991 



C 

C 

C Values for the gas constant, molar masses, calorific values 

C and the density of air conform with ISO 6976 (1995) . 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE SGERG (X2 , X3 , HS , RM, X5 , P,TC,Z,D) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

IF( P.LT. 0.0 .OR. P .GT.12 0.0) STOP ' PRESSURE OUT OF RANGE' 

IF(TC .LT. -23 . .OR. TC . GT .65.0) STOP ' TEMPERATURE OUT OF RANGE' 

CALL SGERGl ( P , TC , X2 , X3 , X5 , HS , RM, Z , D) 

RETURN 

END 
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RM.GT. 0.90) STOP 
X3.GT. 0.30) STOP 
HS.GT.4 8.0 ) STOP 
4 5*X5) .GT.RM)STOP 



'REL. MASS OUT OF RANGE' 
'C02 OUT OF RANGE' 
'CALOR. VALUE OUT OF RANGE' 
'CONFLICTING INPUT' 



(2 ****************************************************************^t*** 

SUBROUTINE SGERGl { P , TC , Q2 , Q3 , Q5 , QM , RM , Z , D) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
COMMON /RBLOK/ AM0L,HS 

COMMON /XBLOK/ XI , X2 , X3 , Xll , X12 , X13 , X22 , X23 , X33 
> ,X5,X7,X15,X17,X25,X55,X77 

COMMON /MBLOK/ GMIRO , GMIRI , GM2 , GM3 , GM5 , GM7 , FA, FB, RL, TO , H5 , H7 , R 
HS = QM 
X3 = Q3 
X5 = Q5 
C IF(RM.LT. 0.5 5. OR 
C IF(X3.LT. 0.0 .OR 
C IF(HS.LT.20.0 .OR 
IF( (0.55+0.97*X3-0 
SM = RM*RL 
X7 = X5*0.0964D0 
X3 3 = X3*X3 
X55 - X5*X5 
X7 7 = X7*X7 
BEFF= -0.065D0 
H = 1000. ODO 
AMOL= l.ODO/ (FA+BEFF) 
K = 
KK = 
1 CALL SMBER(H, SMTl) 

IF(ABS(SM-SMT1) . GT . l.D-5) THEN 

CALL SMBER(H+1.0D0,SMT2) 

DH= (SM-SMTl) / {SMT2-SMT1) 

H = H+DH 

KK = KK+1 

IF(KK.GT.20)STOP ' NO CONVERGENCY 

GO TO 1 
END IF 
Xll = X1*X1 
X12 = X1*X2 
X13 = X1*X3 
X22 = X2*X2 
X2 3 = X2*X3 
X2 5 = X2*X5 
X15 = X1*X5 
X17 = X1*X7 
CALL B11BER(T0,H,B11) 
CALL BBER{T0,B11,BEFF) 
AMOL= 1. ODO/ (FA+BEFF) 

HSBER = Xl*H*AMOL+(X5*H5+X7*H7) *AMOL 
IF{ABS(HS-HSBER) .GT.l.OD-4) THEN 

K = K+1 



#1 



IFfK.GT. 


.20) STOP 


'NO 


CONVERGENCY 


#2' 




GO TO 1 
















END IF 
















IF{X2 .LT. -0, 


.01 . 


OR. X2 


.GT 


.0.5) 


STOP 


'CALC 


. N2 


IF(X2+X3 .GT 


.0.5) 


STOP 




■N2 + 


C02 OUT 


OF RANGE' 


IF(0.55+.4*X2+0. 


,97*X3- 


0.4 


5*X5) 


-GT.RM) 










STOP 




CONFLICTING RESULT 


FOR 


Q2 = X2 
















T = TC+TO 
















CALL B11BER(T,H, 


BID 












CALL BBER(T, 


,B11, 


B) 













OUT OF RANGE' 



N2' 



CALL CBER(T,H,C) 

CALL ITER(P,T,B,C,V,Z) 

D = l.ODO/V 

END 

SUBROUTINE SMBER(H,SM) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /RBLOK/ AM0L,HS 

COMMON /XBLOK/ XI , X2 , X3 , Xll , X12 , X13 , X22 , X23 , X33 

,X5,X7,X15,X17,X2 5,X5 5,X77 
COMMON /MBLOK/ GMIRO , GMIRI , GM2 , GM3 , GM5 , GM7 , FA, FB, RL, TO , H5 , H7 , R 
GM1= GM1R0+GM1R1*H 
XI = (HS-(X5*H5+X7*H7) *AMOL) /H/AMOL 
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X2 = 1.0D0-X1-X3-X5-X7 

SM = (X1*GM1+X2*GM2+X3*GM3+X5*GM5+X7*GM7) *AMOL 

END 

***-k********-ki<i<*****-k******************-**:******* ******************** 

SUBROUTINE BllBER (T, H, Bll ) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /BBLOK/ BR11H0(3), BRllHl(3), BR11H2(3), BR22(3), BR23{3), 
+ BR33(3), BR15(3), BR17(3), BR55(3), BR77(3), B25 

T2=T*T 

B11=BR11H0(1) + BR11H0{2)*T + BR11H0(3)*T2 
+ +(BR11H1(1) + BR11H1{2)*T + BRllHl (3 ) *T2 ) *H 
+ +(BR11H2(1) + BR11H2{2)*T + BR11H2 ( 3 ) *T2 ) *H*H 

END 

SUBROUTINE BBER (T, Bll , BEFF) 
IMPLICIT DOUBLE PRECISION (A-H/O-Z) 

COMMON /BBLOK/ BR11H0(3), BR11H1(3), BR11H2(3), BR22{3), BR23(3), 
+ BR33(3), BR15(3), BR17(3), BR55(3), BR77(3), B25 

COMMON/ZETA/ Z12 , Z13 , Y12 , Y13 , Y123 , Y115 
COMMON /XBLOK/ XI , X2 , X3 , Xll , X12 , X13 , X22 , X23 , X33 

> ,X5,X7,X15,X17,X25,X55,X77 
T2=T*T 

B22=BR22{1) + BR22(2)*T + BR22(3)*T2 

B23=BR23(1) + BR23(2)*T + BR23(3)*T2 

B3 3=BR33(1) + BR33(2)*T + BR33(3)*T2 

B15=BR15(1) + BR15(2)*T + BR15(3)*T2 

B55=BR55(1) + BR55(2)*T + BR55(3)*T2 

B17=BR17(1) + BR17(2)*T + BR17(3)*T2 

B77=BR77(1) + BR77(2)*T + BR77(3)*T2 

BA13= B11*B33 

IF (BA13 .LT. 0.0) STOP ' NO SOLUTION' 

ZZZ=Z12+ {32 .ODO-T) **2*1.87 5D-5 

BEFF = X11*B11 + X12*ZZZ* (B11+B22) + 2 . 0D0*X13 *Z13 *DSQRT {BA13 ) 

> +X22*B22 + 2.0D0*X23*B23 + X33*B33 + X55*B55 

> +2 .0D0*X15*B15 + 2.0D0*X25*B25 + 2 . 0D0*X17*B17 + X77*B77 
END 
*-******************************************************************* 

SUBROUTINE CBER (T, H, CEFF) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /CBLOK/ CR111H0(3), CR111H1(3), CR111H2(3), CR222(3), 
+ CR223(3), CR233(3), CR333(3), CR555(3), CR117(3) 

COMMON /ZETA/ Z12 , Z13 , Y12 , Y13 , Y123 , Y115 
COMMON /XBLOK/ XI , X2 , X3 , Xll , X12 , X13 , X22 , X23 , X33 

> ,X5,X7,X15,X17,X25,X55,X77 
T2=T*T 

Clll^CRlllHO (1) + CR111H0(2)*T + CRlllHO ( 3 ) *T2 
+ +(CR111H1(1) + CR111H1(2)*T + CRlllHl ( 3 ) *T2 ) *H 
+ +{CR111H2(1) + CR111H2(2)*T + CR111H2 ( 3 ) *T2 ) *H*H 

C222 = CR222(1) + CR222(2)*T + CR222(3)*T2 

C223 = CR223(1) + CR223(2)*T + CR223(3)*T2 

C233 = CR233(1) + CR233(2)*T + CR233(3)*T2 

C333 = CR333(1) + CR333(2)*T + CR333(3)*T2 

C555 = CR555(1) + CR555(2)*T + CR555(3)*T2 

C117 = CR117(1) + CR117(2)*T + CR117(3)*T2 

CA112=C111*C111*C222 

CA113=C111*C111*C3 3 3 

CA12 2=C111*C2 2 2*C22 2 

CA12 3=C111*C22 2*C33 3 

CA13 3=C111*C333*C3 3 3 

CA115=C111*C111*C555 

IF {CA112.LT.0. .OR. 
+ .OR. CA123.LT.0.0 
+ STOP ' NO SOLUTION' 

D3REP=1.0D0/3.0D0 

CEFF=X1*X11*C111 +3.0D0'' 
+ *0.0013D0) 

+ + 3.0D0*X11*X3 * (CA113) **D3REP *Y13 
+ + 3.0D0*X1*X15 * (CA115) **D3REP *Y115 
+ + 3.0D0*X1*X22 * (CA122) **D3REP * ( Y12+ (T-270 . ODO ) *0. G013D0) 



CA113 , 


. LT . . 


.OR. 


CA122,LT.0 


.0 




.OR. 


CA13 3 .LT 


.0.0 


.OR 


. CA115. 


LT. 


0.0 


X11*X2* (CA112) 


**D3REP* 


(Y12+(T- 


270. 


ODO) 
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+6 . 0D0*X1*X2*X3* (CA123) **D3REP *Y123 

+3 . 0D0*X1*X33 * {CA133) **D3REP *Y13 
h +X22*X2*C222 + 3 . 0D0*X22 *X3 *C223 + 3 . 0D0*X2 *X33 *C233 
h +X3*X33*C333 + X5*X55*C555 + 3 . 0D0*Xll*X7*C117 

RETURN 
END 

SUBROUTINE ITER ( P, T, B , C, V, Z ) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /MBLOK/ GMlRO , GMlRl , GM2 , GM3 , GM5 , GM7 , FA, FB, RL, TO , H5 , H7 , R 



RT 


= R*T 








RTP= RT/P 








V 


= RTP+B 








KK 


= 








V 


= RTP*(1. 


0D0+B/V+C/V**2) 




KK 


= KK+1 








IF 


( KK . GT . 


20) STOP 


' NO 


COl 


Z 


= 1 . 0D0+B/V+C/V**2 






PA 


:--: RT/V*Z 








IF 


(ABS (PA-P) 


.GE. l.D- 


-5)G0 


TO 


RETURN 









END 



TirTt:A-***-t***-Ar** + ****-Jr****************Tlr********** + *******************: 



BLOCK DATA 

IMPLICIT DOUBLE PRECISION 

COMMON /BBLOK/ BR11H0(3), 



(A-H,0-Z) 

BR11H1(3), BR11H2(3), BR22(3), BR23(3) 
BR33(3), BR15(3), BR17(3), BR55(3), BR77{3), B25 
COMMON /CBLOK/ CR111H0(3), CRlllHl(3), CR111H2(3), CR222(3), 



+ 




CR223(3) , 


, CR233(3), CR333(3), CR555(3), CR117(3 


COMMON /ZETA/ 


Z12,Z13,Y12,Y13,YI2 3, Y115 


COMMON /MBLOK/ GMIRO , GMIRI , GM2 , GM3 , GM5 , 


, GM7 , FA , FB , RL , TO , H5 , H7 , R 


DATA 


BRllHO/- 


0.425468D0, 


0.286500D-2, ■ 


■■.462073D-5 /, 


+ 


BRllHl/0 


.877118D-3, 


-.5552810^5, ( 


3.881510D-8 /, 


-(- 


BR11H2/- 


.824747D-5, 


0.431436D-8, - 


- .608319D-1I/, 


+ 


BR22 / - 


. 144600D0 , 


0.740910D-3, - 


-.911950D-6 /, 


+ 


BR23 / - 


.339693D0 , 


0.161I76D-2, - 


-.204429D-5 /, 


+ 


BR3 3 / - 


.868340D0 , 


0.403760D-2, ■ 


-.516570D-5 /, 


+ 


BR15 / - 


.521280D-1, 


0.271570D-3, - 


-.25 ■ D-6 /, 


+ 


BR17 / - 


.687290D-1, 


-.239381D-5, ( 


3.518195D-5 /, 


+ 


BR55 / - 


.110596D-2, 


0.813385D-4, - 


-.987220D-7 /, 


-*- 


BR77 / - 


.130820D0 , 


0.602540D-3, - 


-.644300D-6 /, 


+ 


B25 / 


. 012D0 / 






DATA 


CRlllHO/ 


-.302488D0 , 0.19586ID-2, 


, -.316302D-5 /, 


+ 


CRllIHl/ 


0.646422D- 


-3, -.422875D-5, 


, 0.588157D-8 /, 


+ 


CR111H2/ 


-.332805D- 


-6, 0.223160D-8, 


, - .367713D-11/, 


+ 


CR222 / 


0.784980D- 


'2, -,398950D"4, 


, 0.611870D-7 /, 


+ 


CR223 / 


0.552066D- 


^2, -.168609D~4, 


, 0.157169D-7 /, 


+ 


CR233 / 


0.358783D- 


-2, 0.805674D-5, 


, -.325798D-7 /, 


+ 


CR3 3 3 / 


.2051300^ 


-2, 0.348880D-4, 


, -.837030D-7 /, 


+ 


CR5 55 / 


0.104711D- 


■2, -.364887D-5, 


.467095D-8 /, 


+ 


CR117 / 


0.736748D- 


-2, -.276578D--4, 


.343051D-7 / 


DATA 


Z12 / 


0.72D0 /, 


Z13 / -.865D0/, 




-+- 


Y12 / 
Y115/ 


0.92D0 /, 
1.2DG/ 


Y13 / 0.92D0/, 


, Y123 /I.IODO/, 


DATA 


GM1R0/-2 


.709328D0/, 


GMIRI/ .02I062199DO/, 


+ 


GM2 / ; 


28.01-35D0/, 


GM3 / 44.010D0/, 


+ 


GM5 / 


2.0159D0/, 


GM7 / 28. ( 


DIODO/, 


+ 


FA /22 


.414097D0/, 


FB / 22.71081ID0/, 


+ 


RL /I 


.292923D0/, 


TO / 273 


.15D0/, 


+ 


H5 / 


285.83D0/, 


H7 / 282 


.98D0/ 


+ 


R /0.( 


3831451D0/ 







END 
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